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Abstract 


A  closed  form  solution  for  the  flow  of  resin  in  the  vacuum-assisted  resin  transfer 
molding  (VARTM)  process  is  used  extensively  for  affordable  manufacturing  of  large 
composite  structures.  During  VARTM  processing,  a  highly  permeable  distribution 
medium  is  incorporated  into  the  preform  as  a  surface  layer.  During  infusion,  the  resin 
flows  preferentially  across  the  surface,  simultaneously  through  the  preform,  to  a  complex 
flow  front.  The  analytical  solution  presented  here  provides  insight  into  the  scaling  laws 
governing  fill  times  and  resin  inlet  placement  as  a  function  of  the  properties  of  the 
preform,  distribution  media,  and  resin.  The  formulation  assumes  that  the  flow  is  fully 
developed  and  is  divided  into  two  areas:  (1)  a  saturated  region  with  no  crossflow,  and  (2) 
a  flow  front  region,  which  moves  with  a  uniform  velocity,  where  the  resin  is  infiltrating 
into  the  preform  from  the  distribution  medium.  The  law  of  conservation  of  mass  and 
Darcy’s  Law  for  flow  through  porous  media  are  applied  in  each  region.  The  resulting 
equations  are  nondimensionalized  and  are  solved  to  yield  the  flow  front  shape  and  the 
development  of  the  saturated  region.  It  is  found  that  the  flow  front  is  parabolic  in  shape, 
and  the  length  of  the  saturated  region  is  proportional  to  the  square  root  of  the  time 
elapsed.  The  obtained  results  are  compared  to  data  from  full-scale  simulation  and  show 
good  agreement.  The  solution  allows  greater  insight  into  the  physics  process,  enables 
parametric  and  optimization  studies,  and  can  reduce  the  computational  cost  of  full-scale, 
three-dimensional  (3-D)  simulations. 
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1.  Introduction 


The  vacuum-assisted  resin  transfer  molding  (VARTM)  process  offers  numerous  cost 
advantages  over  traditional  RTM  via  lower  tooling  costs,  room  temperature  processing,  and 
scalability  to  large  structures.  Recent  advanced  technology  demonstrators  such  as  the  Advanced 
Enclosed  Mast  Sensor  (AEM/S)  System  and  the  composite  advanced  vehicle  (CAV)  have  shown 
the  potential  of  VARTM  technology  for  the  low-cost  fabrication  of  large-scale  structures 
requiring  thick-section  construction  and  hybrid  multifunctional  integral  armor.  The  VARTM 
process  is  also  used  extensively  in  commercial  applications  such  as  bridge  decks,  rail  cars,  and 
yachts  (Figure  1). 


Figure  1.  Examples  of  the  Broad  Application  Potential  for  VARTM  Processes  Including 
Shipping,  Infrastructure,  Land  Combat  Vehicle  Armor,  and  Repair. 


1 


VARTM  is  a  composites  manufacturing  process  that  involves  the  layup  and  vacuum  bagging 
of  dry  reinforcing  fibers  in  fabric,  tape,  or  bulk  form  as  a  preform  in  a  one-sided  open  mold  and 
impregnating  the  preform  with  liquid  resin  using  negative  pressure  (i.e.,  a  vacuum)  followed  by 
cure  and  demolding.  The  advantages  of  the  VARTM  process  over  the  RTM  process  are 
scalability  and  affordability  for  the  fabrication  of  large  composite  structures.  Large  parts  can  be 
infused  rapidly  using  vinyl  esters,  phenolics,  and  epoxies  at  room  temperature  under  vacuum 
pressure  only.  Consequently,  tooling  costs  and  investments  are  substantially  reduced.  VARTM 
is  a  completely  closed  system  that  traps  volatile  organic  compounds  (VOCs),  reduces  the  need 
for  solvents,  and  results  in  less  scrap  than  other  processes. 

The  present  study  focuses  on  Seemann’s  Resin  Infusion  Molding  Process  (SCRIMP) 
(Seemann  1990).  In  this  VARTM  process,  a  highly  permeable  distribution  medium  is 
incorporated  into  preform  as  a  surface  layer.  During  infusion,  the  resin  flows  preferentially 
across  the  surface  and  simultaneously  through  the  preform  thickness  enabling  large  parts  to  be 
fabricated  solely  under  vacuum  pressure.  The  layup  of  the  materials  in  the  process  is  shown  in 
Figure  2. 

In  very  large  composite  structures,  multiple  inlet  gates  are  required  to.  ensure  complete  wet- 
out  of  the  part  prior  to  gelation  of  the  resin.  Selection  of  distribution  media,  performs,  and  gate 
and  vent  locations  are  based  on  past  experience  for  similar  applications.  New  applications  in 
which  part  thickness,  resin,  or  preform  characteristics  change  require  costly  trial  and  error 
process  development.  Hence,  a  fundamental  understanding  of  the  process  physics  and  associated 
models  represent  a  significant  contribution  to  the  science  base  for  VARTM. 

Modeling  and  predicting  the  flow  during  the  impregnation  process  provides  insight  into  the 
process  physics  and  highlights  potential  problems  before  production.  In  addition,  flow 
prediction  enables  optimization  of  the  design  variables  affecting  the  process,  such  as  the  distance 
between  resin  inlets  (in  the  case  of  multiple  lines  and  thickness  of  the  diffusion  layer),  and 
provides  rules  of  thumb  for  scaling  of  the  prototypes.  Thus,  a  fundamental  understanding  of  the 
underlying  science  will  help  develop  models  to  reduce  costs,  aid  in  selection  of  design 
parameters,  and  improve  quality. 
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The  flow  of  resin  through  porous  media  such  as  fiber  preforms  and  resin  distribution  media  is 
governed  by  Darcy’s  Law: 


«=— -VP,  (1) 

where  u  is  the  Darcy’s  velocity  (defined  as  the  total  flow  rate  per  total  flow  front  area),  K  is  the 
permeability  tensor  (which  characterizes  the  ease  of  flow  through  the  fiber  perform),  and  /u  is 
the  viscosity  of  the  resin.  This,  when  coupled  with  the  continuity  equation  for  incompressible 
flow,  gives  the  Laplace  equation  for  the  fluid  pressure  field  inside  a  region  permeated  by  the 
fluid: 


V- 


=  0 . 


(2) 


This  equation  can  be  discretized  using  finite  element  methods  which  then  form  the  basis  for 
simulation  of  mold-filling  during  the  resin  infusion  process  (Bruschke  and  Advani  1990,  1991a; 
Liu  et  al.  1996;  Mohan  et  al.  1999). 

The  flow  simulations  can  be  either  two-dimensional  (2-D)  or  three-dimensional  (3-D).  In 
2-D  flow  modeling  (Bruschke  and  Advani  1991b;  Trochu  et  al.  1994;  Lee  et  al.  1994),  the  flow 
of  resin  through  the  thickness  is  considered  uniform,  and  the  finite  element  discretization  is 
applied  along  the  other  two  directions  as  with  liquid  injection  molding  simulation  (LIMS),  which 
is  based  on  the  finite  element/control  volume  approach.  In  2-D  simulations,  only  the  in-plane 
permeabilities  are  supplied  (Simacek  et  al.  1998).  In  3-D  simulations,  the  pressure  and  flow  in 
all  three  directions  is  solved,  and  a  3-D  permeability  sensor  is  supplied  as  input,  as  in  the  resin 
infusion  process  simulation  (REPS),  which  is  based  on  finite  element  methods  without  the  use  of 
the  control  volume  approach  (Gallez  and  Advani  1996).  Usually,  the  geometry,  the  material 
parameters,  and  the  position  of  resin  inlets  and  outlets  are  specified  before  the  filling  simulation 
is  carried  out.  Simulation  codes  are  used  to  track  flow  fronts  and  estimate  the  fill  times. 
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Parametric  studies  then  can  be  conducted  with  simulations  to  design  the  mold  and  the  process 
parameters. 

Closed  form  analytical  solutions  have  also  been  derived  for  the  resin  flow  under  simplifying 
assumptions  and  for  simple  geometries.  These  solutions  explain  the  role  of  various  process 
variables  and  their  interactions  during  processing.  Indeed,  a  closed  form  solution  of  the  resin 
flow  during  the  VARTM  process  not  only  enables  parametric  studies,  optimization,  and 
reduction  of  computational  expenses  of  full-scale  simulations,  but  also  offers  insight  on  the 
scaleup  of  the  process  and  material  parameters  for  large  structures. 

In  earlier  work,  Tari  et  al.  (1998)  derived  a  closed  form  model  for  vacuum-bag  RTM  under 
several  simplifying  assumptions.  They  assumed  that  the  velocity  of  resin  in  the  fiber  preform  is 
negligible  and  that  the  region  behind  the  flow  front  is  uniformly  saturated.  In  the  present  work, 
these  assumptions  were  not  made,  so  the  velocity  of  the  resin,  as  well  as  the  shape  of  the  flow 
front  through  the  thickness  of  the  fiber  perform,  are  accurately  captured.  This  is  important  for 
scaling  purposes. 

A  closed  form  solution  for  the  flow  of  resin  in  the  VARTM  process  is  presented  here.  The 
layup  is  modeled  as  the  distribution  layer  (high  permeability  material)  and  the  structural  layer 
(preform  material).  It  is  assumed  that  the  flow  is  well  developed  and  can  be  divided  into  a 
saturated  region  with  no  crossflow  and  a  flow  front  region  in  which  the  resin  infiltrates  the 
preform  from  the  distribution  medium.  The  flow  front  region  is  assumed  to  be  fully  developed 
with  a  uniform  velocity.  The  law  of  conservation  of  mass  and  Darcy’s  Law  for  flow  through 
porous  media  are  applied  in  each  region.  The  resulting  equations  are  nondimensionalized  and 
are  solved  to  yield  the  flow  front  shape  and  the  development  of  the  saturated  region. 

2.  Problem  Statement 

As  illustrated  in  Figure  3,  the  layup  of  materials  is  modeled  as  two  layers  of  permeable 
materials.  The  distribution  layer  is  much  thinner  than  the  structural  layer,  h{  «  h2 ,  where  hx  and 
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PQ :  Injection  Pressure 


Note:  d  -  region  with  transverse  flow  =  flow  front  length. 

D  -  region  without  transverse  flow  =  length  behind  the  flow  front  region. 

UF  -  flow  front  velocity. 

fi  -  viscosity  of  resin. 

Figure  3.  Two-Layer  Model  of  Resin  Flow  in  the  VARTM  Process. 

are  the  respective  thicknesses  of  the  two  layers.  The  flow  front  in  the  distribution  layer  is 
considered  uniform.  The  permeability  of  the  distribution  layer  is  Klxx  along  the  flow  direction, 
and  the  permeabilities  of  the  structural  layer  are  K2xx  and  K2yy  in  the  x  and  y  directions, 
respectively.  The  constant  inlet  injection  pressure  (atmospheric  pressure)  is  P0 ,  and  the  resin 
viscosity  is  ju . 

In  the  saturated  region,  the  flow  is  one-dimensional  (1-D)  with  Darcy’s  velocities  Ux  and  U2 
in  layers  1  and  2,  respectively.  The  length  of  this  saturated  region  is  D,  and  the  pressure  at  its 
boundary  with  the  second  flow  region  is  assumed  to  beP^ . 

The  second  region,  illustrated  in  Figure  4,  is  the  flow  front  region  where  there  is  transverse 
flow  from  the  distribution  layer  to  the  structural  layer.  The  flow  front  region  of  length  d  is 
assumed  to  maintain  its  shape,  given  by  hF ,  and  advances  with  a  uniform  horizontal  velocity  of 
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Flow  Front  Region 


> 


Figure  4.  Schematic  of  Resin  Flow  in  the  Flow  Front  Region  for  the  Two-Layer  VARTM 
Model. 

UF .  This  is  the  observed  velocity  of  the  resin  and  not  the  Darcy’s  velocity.  The  transverse 
velocity  of  resin  infiltration  from  the  distribution  layer  into  the  structural  layer  is  ul2y .  The 
horizontal  velocity  in  the  flow  front  region  in  the  distribution  layer  ulx  is  with  boundary 
condition  ulx(x  =  D  +  d)  =  0lUF. 

Since  the  resin  is  an  incompressible  fluid,  using  the  continuity  equation  and  Darcy’s  Law  in 
the  structural  layer,  the  governing  equation  for  the  pressure  distribution  is 
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(3) 


K 


d2P 


+  K, 


d  P 


2xx  8x 2  2yy  dy 2 


=  0. 


Consider  the  following  nondimensional  variables: 


D*  P  *  x  .  y 

P  =—  =— =  — > 

pc  xc  y c 


(4) 


where  xc  and  yc  are  the  characteristic  length  scales  in  the  longitudinal  and  thickness  directions, 

respectively,  and  Pc  is  a  scaling  parameter  for  the  pressure.  Introducing  the  dimensionless 

variables,  the  governing  equation  can  be  recast  in  dimensionless  form  as  follows  (Pillai  and 
Advani  1998a,  1998b): 


rr  v2  a  *2  a.  *2  u' 

JXlyyXc  OX  dy 


(5) 


Since  the  resin  distribution  media  was  used  in  the  process  to  enable  the  rapid  and  uniform 
distribution  of  the  resin  across  the  mold  surface,  it  can  be  assumed  that  within  the  flow  front 
region  in  the  structural  layer,  the  major  portion  of  the  resin  flow  is  from  the  distribution  layer 
into  the  structural  layer.  Hence,  the  flow  rate  in  the  y  direction  must  be  more  significant  than 
that  in  the  x  direction  (i.e.,  Qy  »  Qx).  The  xc  in  the  flow  front  region  is  d,  while  yc  is  h2 .  If 
u  and  v  are  the  average  Darcy’s  velocities  in  the  two  directions,  then 


Qy  =Vd, 

and  (6) 

Qx  =uh2 . 

Considering  Darcy’s  equation  for  the  velocities,  the  following  scaling  argument  can  be  made: 
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u  ~ 


(7) 


K**SP 
fj.  dx 


K^Pc 

V  d  ’ 


and 


^2yy  dP  .  ^  ^  ^~2w 

M  fy  P  K 


Since  Q  »  Qx ,  from  equations  6  to  8,  it  can  be  determined  that: 


Kiyydl 


«1. 


(8) 


(9) 


This  allows  for  the  neglect  of  the  x-term  in  pressure  equation  5  in  the  flow  front  region,  which 
leads  to  the  following  result: 


d2P 

dy 2 


=  0 


dP 

dy 


=  />(*)> 


and 


(10) 


V  = 


K2yy  dP 


li  dy 


^2  yy 


/f  W- 


At  the  top  of  the  structural  layer,  where  the  flow  is  always  from  the  distribution  layer,  the 
boundary  condition  is  defined  as  v  1^=  uUy{x).  Hence,  in  the  flow  front  region. 


v(x)  =  «i2>,(x), 
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where  uny(x)  is  the  velocity  of  the  resin  flow  from  the  distribution  to  the  structural  layer  in  the 

flow  front  region.  In  the  saturated  region,  the  length  scale  of  the  flow  is  D  in  the  x  direction  and 
h2  in  the  y  direction,  where  D»h2.  Since  D»d,  from  equation  9, 


K2ryd‘  K^D1 


«1. 


(11) 


So,  the  first  term  in  equation  5  can  be  neglected,  and 


d2P 
dy 2 


=  0: 


dP 

dy 


=  fs(x). 


(12) 


Thus,  the  velocity  in  they  direction  in  the  saturated  region  is 


K2yy  dP 
fi  dy 


^2yy 


/,(*)• 


(13) 


At  the  bottom  of  the  structural  layer  in  the  saturated  region,  the  resin  is  in  contact  with  the 
surface  of  the  mold,  which  is  impermeable.  Hence,  the  no-penetration  boundary  condition  was 
applied  (i.e.,  v  =  0  @y  =h2)  in  the  saturated  region.  Therefore, 

dP 

—  =  fs  (v)  =  0  in  the  saturated  region.  (14) 

dy 


As  a  result,  v  =  0  in  the  saturated  region  everywhere  in  the  structural  layer.  Since  v  =  0,  the 
second  term  in  pressure  equation  3  becomes  zero,  and 


as 


dy 


So,  it  can  be  inferred  that  g(y)  is  constant,  and 


u  =  — 


p.  dx 


(15) 


is  constant  in  the  saturated  region.  Hence,  from  boundary  conditions  of  P  =  P0  @x  =  0  and 
P  =  PD@x  =  D, 


Qp  p  _  p 

—  =  — - -  in  the  saturated  region. 

dx  D 


(16) 


3.  Analytical  Solution 


Considering  the  element  fluid  volumes  shown  in  Figure  5  and  invoking  the  mass  balance. 


—  dh 

unrdx  +  U2(-dhF)=UF02(-dhF)^  uUy  =-J-(02Uf  -U2). 


(17) 


A  lumped  mass  balance  in  the  distribution  layer  in  the  flow  front  region  gives 


-  dulxh}  =  un dx  =>  un  =  -A, 


du 


lx 


dx 


(18) 


11 


dx 


< - ► 

- ► 

Flow  Front  Region 


Figure  5.  Illustration  for  Resin  Mass  Balance  in  the  Flow  Front  Region  in  the  Two-Layer 
Model. 

From  equations  17  and  18,  the  result  is 


(19) 


Applying  Darcy’s  Law  in  the  y  direction  in  the  flow  front  region  and  using  equations  17  and 


18, 


UX2y  = 


ha-Oz  0 

n  hr 


and 


*ar  = -K ^ -u2)=^-5- 

dx  dx  (i  hF 


(20) 


12 


This  leads  to  the  following  equation  for  P,(jc’)  ,  the  pressure  field  in  the  distribution  medium  in 
the  flow  front  region: 


(21) 


A  similar  mass  balance  in  the  flow  front  region,  including  both  the  structural  and  the  distribution 
layers,  yields 

Ufa  +  U2h2  =UF(0lhl+02h2).  (22) 

Applying  Darcy’s  Law  in  the  saturated  region  with  1-D  flow, 

TT  —  ^]xx  —  ~  1-JL 

1  fi  D  ’ 

and  (23) 


U2 


^2xx  Po  £d 
fi  D 


K, 


\xx 


K 


2  xx 


In  combination  with  equation  22, 


Ui 

UF 


{0^+®2h2) 


K 


\ ' 


K. 


—  h,  +K 


2xx 


(24) 


The  previous  set  of  equations  can  be  nondimensionalized  using  the  following  nondimensional 
variables: 
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tr.SLx.JL 

UF 


x-D 

~hT’ 


(25) 


and 


This  gives  the  following  system  of  equations: 


.  -IdP, 


Ulx  = 


(i  dx 


*  _  {<f>X  +<p2) 
h 


u>=f 


K 

yIX\xx 


■  +  1 


duu  _  (<t>2  -  U*2)dh ; 

”  a; 


(26) 


(27) 


(28) 


and 


dx 


K 


2yy 


(29) 


Here,  fi'  = 


is  obtained  from  the  nondimensional  Darcy’s  Law.  The  boundary  conditions 
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on  u]x  are  u[z (o)  =  Ux,u’x id' )  =  ^,  •  The  boundary  conditions  on  /*;(0)  =  U;(<f  )=0.  The 


pressure  boundary  conditions  are 


p;(o)=p;,p;(d-)=o. 


Integrating  equation  29  and  applying  the  boundary  condition  w*x(o)  =  U’ ,u‘x(d'  )=  <t>x , 


<u 

«i 


+  0, 


(30) 


In  combination  with  equation  26, 


Combining  equations  29  and  31,  the  following  nonlinear  ODE  with  boundary  conditions 
tiF (o)  =  1,  hF  {d * )  =  0  result  in: 


d 

dx' 


dhp 

dx  ’  j 


02-u; 


(32) 


This  ODE  is  of  the  form  given  by  (y2 -  ay  +  b.  It  can  be  solved  using  the  substitution 
dy  dp 

p  =  — ,y  =  p —  (Murphy  1960).  This  usually  yields  a  solution  in  a  quadratic  form. 
dx  dy 

The  quadratic  form  h’F  (x* )  =  ax'2  +  /3x*  +  7  can  be  substituted  into  the  above  equation  to  find  the 
solution  using  the  first  boundary  condition  hF  (0)  =  1  and  matching  the  coefficients  of  the  powers 
of  x*  on  either  side.  Then  the  second  boundary  condition,  hF(d')  =  0,  can  be  used  to  determine 

d' .  This  gives  a  quadratic  equation  for  d*  having  two  roots.  Both  roots  are  positive,  but  if  the 
larger  one  were  chosen,  then  the  flow  front  profile  would  be  physically  impossible  (Figure  6). 
Hence,  the  smaller  root  gives  the  following  unique  solution: 
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a  = 


+ 


OxK 


lyy 


02-U2 


(33) 


and 


7=1. 


Thus, 


h'F(x)  = 


K. 


lyy 


2  K 


6h, 


x 


lyy 


3  h’ 


+ 


02-u; 


X  +1 


(34) 


and  since  hF(d')  =  0,d *  can  be  determined  as 


3h[  f  10,  ~  1 

^{i02-u;+3h; 


g 

0,  -u: 


(35) 


Substituting  the  form 


tiF(x )  =  ax'2  +  j8x*  +  7  in 


equation  29,  the  pressure  can  be 


determined  as  follows: 


(36) 


Applying  the  boundary  condition  P’  (0)  =  P'D  gives 
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*  * 
x  -a 


Nonphysical  Region 


Flow  Front  Shape:  hF(x  ) 


Figure  6.  Two  Mathematical  Roots  of  hp*(x*)  =  0.  Note  That  Only  the  Smaller  Root  Is 
Physically  Possible. 


.  (*,-!/;)»•  (<p,-c/;V  l2x;„  w 

D  rjr*  P.  -rr*  nf  ~ 


3 h;  02  -u; 


Using  the  nondimensional  form  of  Darcy’s  Equation  in  the  saturated  region, 


U'2  =^-LA  =>£/*£>*  = 


. 02~U 2*  \lK^~0^2yy 


H  D 


K2yy  V  3A,  <Z>2-£/2 


To  find  C4  substituted  with  fi  is 


jjh2Uf_  =  _]__ 
kXxxp0  ro‘  +A 
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where 


r  = 


and 


A  - 


&-u\ _) 


K 


2yy 


'll 


2K 


2yy  +  yy 


3  hi  02-U2 


From  equation  38,  the  flow  front  velocity  is 


UF  = 


KuA  1 

I uh2  ID *  +yl 


(39) 


(40) 


but  UF  =  ,  and  D'  =  — ;  hence,  the  previous  equation  becomes 

dt  h2 


dD  K+A  1 

dt  fjth2  ID  +  A  h2 


Solving  the  resulting  differential  equation  for  D(t), 

t-tQ=cXD2-Dl)+C2{D-D,), 


(41) 


(42) 


where  C,  =  — — —  and  C,  =  ^  — 
2KlxxP0  K1XXP0 


The  variable  t0  is  the  time  it  takes  for  the  flow  front 


region  to  become  fully  developed,  while  D0  is  the  entry  length  for  the  development  of  the  flow 
front  region.  Solving  for  D(t), 
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D(t)  = 


(43) 


(Ah2f  +—~-—(t-t0+C]D*+C2D0)-Ah2 

_ •4**ixr‘o _ 

r 

Equations  42  and  43  are  important  for  the  design  of  the  VARTM  process.  The  results 
obtained  are  compared  to  full-scale,  finite-element-based  simulations  using  LIMS  4.0  and  are 
presented  in  the  next  section.  A  parametric  study  is  presented  to  shed  light  on  how  one  can  scale 
the  parameters  in  the  VARTM  process. 

4.  Verification:  Full-Scale  Simulations 


The  results  obtained  can  be  compared  to  finite-element-based  simulations  of  the  filling 
process  in  VARTM.  LIMS  4.0  was  used  for  simulation  of  the  filling  process  for  five  different 
cases,  with  different  values  for  permeabilities  and  fiber  volume  fractions  for  the  distribution  and 
structural  layers,  respectively,  and  length  of  part,  D.  Each  part  was  modeled  using  finite 
elements  and  the  filling  process  simulated  as  a  constant  pressure  injection  at  1  atm  at  one  comer 
of  the  part.  The  viscosity  of  the  resin,  the  thicknesses  of  the  structural  and  distribution  layers, 
and  the  permeability  of  the  distribution  layer  were  held  constant.  The  fill  times  and  the  values  of 
d  were  determined  and  compared  to  those  obtained  from  the  closed  form  solution.  The  results 
are  tabulated  in  Table  1.  The  fill  time  from  the  analytical  solution  was  found  to  be  within  2%  of 


the  value  from  the  full-scale  numerical  simulation,  while  the  value  of  d  was  within  12%  of  the 


value  from  simulation, 
analytical  solution. 


For  all  cases,  the  condition 


« 1 .  was  maintained  for  a  valid 


The  flow  front  history  and  the  pressure  contours  at  the  final  time  are  plotted  for  Case  1  in 
Figure  7.  It  can  be  observed  that  the  flow  front  is  constant  in  shape,  while  the  lines  of  constant 
pressure  in  the  saturated  region  are  equally  spaced  and  vertical  to  the  x  axis,  thus  verifying  the 
assumptions  of  constant  flow  front  shape  and  linear  variation  in  pressure  in  the  saturated  region. 
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Table  1.  Comparison  of  Closed  Form  Solution  With  Results  From  Full-Scale  Numerical 
Solution  of  the  VARTM  Process 


Case 

Parameters 

d 

(cm) 

t-t0 

(s) 

%  Error3 

LIMS  4.0 

Closed  Form 
Solution 

LIMS  4.0 

Closed  Form 
Solution 

d 

(%) 

1 

1 

8.9 

9.8 

44.9 

44.8 

2 

29.2 

30.9 

19.1 

19.4 

5.7 

1.5 

3 

3.0 

3.1 

54.2 

54.5 

2.9 

0.5 

1 

9.2 

10.1 

71.7 

72.1 

10.2 

0.6 

5 

K2xx=  8.8E-7  cm2 
K2yy  =  4.4E-7  cm2 

4>1  =  0.70,  <t>2  =  0.50 
D-Do  =  40.0  cm 

9.0 

10.0 

44.2 

44.5 

11.6 

0.7 

Note:  The  following  variables  were  used  in  all  cases: 
P0  =  1E6  g/cm-s2, 
jLi  =  1  g/cm-s 
hj  =  0.01  cm, 
h2  =  1  cm, 

Klxx  “  IE-3  cm2,  and 
Do~0  cm  for  all  the  cases. 
a  Based  on  results  from  LIMS  4.0. 
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5.  Effect  of  Process  Variables:  A  Parametric  Study 


The  process  variables  that  affect  the  flow  of  resin  are  broadly  classified  into  geometric 
parameters  (such  as  thickness)  and  material  properties  (such  as  permeability  and  porosity  of  the 
two  layers).  These  process  variables  influence  the  time  to  fill  a  mold  of  a  given  length.  A 
parametric  study  of  these  effects  allows  for  better  design  and  analysis  of  the  VARTM 
manufacturing  process.  In  the  present  section,  the  effect  of  a  number  of  process  variables  on  the 
fill  times  and  flow  velocity  was  studied.  The  baseline  values  used  for  the  study  are: 

P0  =  1  atm,  /X  =  1  cp, 
h,  =1.00  cm,  h2  =  0.025  cm, 

Klxx  =  8.8xlO“4cm2,K2xx  =  8.8x  10"7cm2,K2yy  =  1.47x  10-7cm2,  and 
=  0.99,  $2  =0.50. 

In  the  plots  for  each  parameter,  the  flow  front  velocity  (Uf)  and  the  time  taken  (/)  are  plotted 
against  D  for  different  values  of  the  parameter.  Note  that  the  time  axis  is  reversed;  therefore,  the 
lines  for  t  start  from  zero  at  the  top  of  the  graph. 

5.1  Effect  of  Thickness  Ratios.  The  effect  of  the  ratio  of  the  thickness  of  the  distribution 

h 

medium  to  that  of  the  structural  layer,  h)  -  —  ,  is  considered.  Figure  8  plots  the  flow  front 
velocity,  UF ,  vs.  D  for  different  thickness  ratios.  As  A,  increases,  the  flow  front  velocity 
increases  while  the  fill  time  decreases  for  a  given  length  D.  This  is  because  as  the  thickness  of 
the  highly  permeable  distribution  medium  relative  to  that  of  the  structural  layer  increases,  the 
resin  flow  rate  in  the  distribution  medium  increases.  Since  the  diffusion  material  is  used  to 
distribute  resin  in  the  part  and  ensure  mold  filling,  an  increase  in  h*  will  cause  an  increase  in 
flow  front  velocity  and  a  decrease  in  fill  time. 

For  the  case  of  h’  =  0.1 ,  encountered  with  thin  section  composite  parts,  UF  shows  a  slower 
decrease  with  D  than  with  the  other  cases.  With  thick-section  composite  parts  where  h\  0 , 
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D(cm) 


Figure  8.  Flow  Front  Velocity  and  Fill  Time  as  a  Function  of  Length  of  the  Saturated 
Region:  Effect  of  Thickness  Ratios. 

JJF  falls  rapidly.  In  order  to  fill  such  a  part  efficiently,  the  distance  between  the  gates  (D)  has  to 
be  small;  thus,  the  number  of  gates  required  increases.  Therefore,  the  solution  provides  insight 
into  the  scaling  laws  required  for  manufacturing  thick-section  composites  by  VARTM. 

5.2  Effect  of  Permeability.  The  permeability  values  of  importance  are:  K']xx,  the 
permeability  of  the  distribution  layer  in  the  longitudinal  x  direction,  and  K'2xx  and  K’2yy,  the 
permeabilities  of  the  structural  layer  in  the  longitudinal  x  and  thickness  y  directions.  As  seen  in 
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Figures  9-11,  as  the  permeability  values  increase,  the  time  to  fill  decreases,  and  the  flow  front 
velocity  increases.  As  the  permeability  values  increase,  the  resistance  of  the  material  to  the  resin 
flow  decreases.  Hence,  the  net  flow  rates  are  higher,  and  the  time  to  fill  decreases,  while  the 
flow  front  velocity  increases.  Since  the  flow  rate  in  the  distribution  layer  is  higher,  the  effect  of 
K‘Xxx  is  significantly  more  than  that  of  K'2xx  and  K\  .  These  effects  can  also  be  observed  from 
the  plots  shown.  This  has  important  ramifications  on  the  selection  and  the  thickness  of 
distribution  media. 

5.3  Effect  of  Porosity.  The  porosity  of  a  fiber  preform  is  defined  as  the  fraction  of  the  total 
volume  of  the  material  not  occupied  by  the  fibers.  In  composite  manufacturing,  the 
complementary  term,  volume  fraction,  was  more  commonly  used.  The  volume  fraction  is 
defined  as  the  fraction  of  the  fiber  preform  occupied  by  the  fibers  and  is  related  to  the  porosity 
by  the  relation  Vf  =1-0 .  The  porosity  also  affects  the  permeability  of  the  material.  However, 
this  coupling  has  not  been  accounted  for  in  the  present  work. 

The  porosity  values  considered  here  are  <Z>, ,  the  porosity  of  the  distribution  layer,  and  02 , 
the  porosity  of  the  structural  layer.  As  observed  from  Figure  12,  <2>,  did  not  significantly  affect 
the  time  to  fill  and  the  flow  front  velocity.  In  Figure  13,  as  &2  increases,  the  flow-front  velocity 
decreases,  and  the  fill  time  significantly  increases.  This  is  because  the  fraction  of  the  total  part 
volume  occupied  by  the  thin  layer  of  diffusion  material  is  very  low  compared  to  that  occupied  by 
the  fiber  preform  in  the  structural  layer.  Hence,  increasing  the  porosity  of  the  diffusion  material 
does  not  have  a  significant  effect.  Whereas,  if  <P2  increases,  the  volume  of  the  structural  layer 
(which  is  unoccupied  by  the  fiber  perform)  significantly  increases.  This  additional  volume  must 
be  filled  by  the  resin,  thus  requiring  more  time  to  fill  it  and  therefore  slowing  down  the  flow 
front. 
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D(cm) 


Figure  9.  Flow  Front  Velocity  and  Fill  Time  as  a  Function  of  Length  of  the  Saturated 
Region:  Effect  of  Permeability  of  Distribution  Medium. 


6.  Conclusions 


A  closed  form  solution  for  flow  of  resin  in  the  VARTM  process  has  been  developed.  This 
process  is  explained  by  a  two-layer  model  comprised  of  a  distribution  layer  and  a  structural 
layer,  containing  fiber  preform.  The  flow  is  divided  into  a  saturated  region  where  there  is  no 
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D(cm) 

Figure  10.  Flow  Front  Velocity  and  Fill  Time  as  a  Function  of  Length  of  the  Saturated 
Region:  Effect  of  In-Plane  Permeability  of  Fiber  Preform. 

crossflow,  and  a  flow  front  region  with  a  steady  shape  and  uniform  velocity  where  the  driving 
flow  emanates  from  the  distribution  layer  to  the  structural  layer.  It  is  assumed  that  the  thickness 
of  the  distribution  layer  is  much  smaller  than  that  of  the  structural  layer,  and  that  the  length  of 
the  flow  front  region  is  much  smaller  than  that  of  the  saturated  region.  It  is  also  assumed  that  the 
crossflow  in  the  flow  front  region  is  much  higher  than  the  flow  from  the  saturated  region  into  the 
flow  front  region.  Darcy’s  Law  regarding  flow  in  porous  media  and  mass  balances  at  different 
sections  was  used  to  formulate  a  system  of  differential  equations,  and  a  closed  form  solution  was 
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Figure  11.  Flow  Front  Velocity  and  Fill  Time  as  a  Function  of  Length  of  the  Saturated 
Region:  Effect  of  Transverse  Permeability  of  Fiber  Preform. 

found.  The  model  predicted  the  shape  and  development  of  the  flow  front  given  the  material 
properties,  the  geometric  parameters,  the  pressure  at  the  inlet,  and  the  viscosity.  The  obtained 
results  were  verified  by  comparing  them  with  full-scale  simulations.  The  parametric  study 
indicated  trends  that  reflect  the  physics  of  the  flow  process  and  identified  the  parameters  that 
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D  (cm) 


Figure  12.  Flow  Front  Velocity  and  Fill  Time  as  a  Function  of  Length  of  the  Saturated 
Region:  Effect  of  Porosity  of  Distribution  Medium. 

significantly  affect  the  filling  process.  Subsquent  solutions  provided  physical  insight  into  the 
manufacturing  process  and  can  be  used  for  scaling,  design,  and  optimization  of  the  VARTM 
process. 
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AMSRL  OP  SD  TP 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 

HQDA 
DAMI  FIT 
NOLAN  BLDG 

WASHINGTON  DC  20310-1025 

DIRECTOR 
DA  OASARDA 
SARD  SO 

103  ARMY  PENTAGON 
WASHINGTON  DC  20310-0103 

DPTY  ASST  SECY  FOR  R&T 
SARD  TT 
THE PENTAGON 
RM  3EA79 

WASHINGTON  DC  20301-7100 

COMMANDER 

US  ARMY  MATERIEL  CMD 

AMXMIINT 

5001  EISENHOWER  AVE 

ALEXANDRIA  VA  22333-0001 


2  COMMANDER 

US  ARMY  ARDEC 
AMSTAARAEWW 
E BAKER 
J  PEARSON 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTAARTD 
C  SPINELLI 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTAARFSE 
PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 

US  ARMY  ARDEC 
AMSTAARCCHA 
W  ANDREWS 
S  MUSALLI 
R  CARR 
M  LUCIANO 
E  LOGSDEN 
T  LOUZEIRO 
PICATINNY  ARSENAL  NJ 
07806-5000 

4  COMMANDER 

US  ARMY  ARDEC 
AMSTAARCC 
G PAYNE 
J  GEHBAUER 
C  BAULIEU 
HOPAT 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  CCH  P 
J  LUTZ 

PICATINNY  ARSENAL  NJ 
07806-5000 
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COPIES 


ORGANIZATION 


1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  FSF  T 
C  LIVECCHLA 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  QAC  T  C 
C  PATEL 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARM 

D  DEMELLA 
F  DIORIO 

PICATINNY  ARSENAL  NJ 
07806-5000 

3  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  FSA 
A  WARNASH 
BMACHAK 

M  CHIEFA 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARFSPG 
M  SCHIKSNIS 
D  CARLUCCI 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTAARFSPA 
P  KISATSKY 
PICATINNY  ARSENAL  NJ 
07806-5000 


9  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  CCH  B 
P  DONADIA 
F  DONLON 
P  VALENTI 
C  KNUTSON 
G  EUSTICE 
S PATEL 
G  WAGNECZ 
RSAYER 
F CHANG 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  CCH  C 
HCHANIN 
S  CHICO 

PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  CCL 
F  PUZYCKI 
R  MCHUGH 
D  CONWAY 
E  JAROSZEWSKI 
R  SCHLENNER 
MCLUNE 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  QAC  T 
D  RIGOGLIOSO 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  SRE 
DYEE 

PICATINNY  ARSENAL  NJ 
07806-5000 
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COPIES  ORGANIZATION 


1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  WET 
TSACHAR 
BLDG  172 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ASF 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  US  ARMY  ARDEC 
INTELLIGENCE  SPECIALIST 
AMSTA  ARWELF 

M  GUERRIERE 
PICATINNY  ARSENAL  NJ 
07806-5000 

11  PMTMAS 

SFAE  GSSC  TMA 
R  MORRIS 
C  KIMKER 
D  GUZOWICZ 
E  KOPACZ 
R  ROESER 
R  DARCY 
R  MCDANOLDS 
LDULISSE 
C  ROLLER 
J  MCGREEN 
B  PATTER 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  PEO  FIELD  ARTILLERY  SYS 
SFAE  FAS  PM 

H  GOLDMAN 
T  MCWILLIAMS 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  TACOM 
PM  ABRAMS 
SFAE  ASM  AB 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 


6  PM  SADARM 

SFAE  GCSS  SD 
COL  B  ELLIS 
M  DEVINE 
R  KOWALSKI 
W  DEMASSI 
J  PRITCHARD 
SHROWNAK 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
PRODUCTION  BASE 
MODERN  ACTY 
AMSMC  PBM  K 
PICATINNY  ARSENAL  NJ 
07806-5000 

3  COMMANDER 

US  ARMY  TACOM 

PM  TACTICAL  VEHICLES 

SFAE  TVL 

SFAETVM 

SFAETVH 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

US  ARMY  TACOM 
PM  BFVS 
SFAE  ASM  BV 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

US  ARMY  TACOM 

PMAFAS 

SFAE  ASM  AF 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 

1  COMMANDER 

US  ARMY  TACOM 
PM  RDT&E 
SFAE  GCSS  W  AB 
J  GODELL 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 
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2 


1 


1 


1 


1 


ORGANIZATION 


NO.  OF 
COPIES 


COMMANDER  14 

US  ARMY  TACOM 

PM  SURV  SYS 

SFAE  ASM  SS 

TDEAN 

SFAE  GCSS  W  GSI M 
D  COCHRAN 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 
US  ARMY  TACOM 
PM  SURVIVABLE  SYSTEMS 
SFAE  GCSS  W  GSI  H 
M  RYZYI 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 
US  ARMY  TACOM 
PM  BFV 

SFAE  GCSS  W  BV 

S  DAVIS  1 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 

US  ARMY  TACOM  1 

PM  LIGHT  TACTICAL  VHCLS 

AMSTA  TR  S 

A  J  MILLS  MS  209 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 

COMMANDER  1 

US  ARMY  TACOM 

PM  GROUND  SYSTEMS 

INTEGRATION 

SFAE  GCSS  W  GSI 

R  LABATILLE 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000  2 

COMMANDER 
US  ARMY  TACOM 
CHIEF  ABRAMS  TESTING 
SFAE  GCSS  W  AB  QT 
T  KRASKIEWICZ 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 


ORGANIZATION 

COMMANDER 
US  ARMY  TACOM 
AMSTA  TRR 
J  CHAPIN 
R  MCCLELLAND 
D  THOMAS 
J BENNETT 
D  HANSEN 
AMSTA  JSK 
S  GOODMAN 
J  FLORENCE 
K  IYER 
J THOMSON 
AMSTA  TRD 
D  OSTBERG 
L  HINOJOSA 
BRAJU 
AMSTA  CS  SF 
H  HUTCHINSON 
F  SCHWARZ 
WARREN  MI  48397-5000 

COMMANDER 
US  ARMY  TACOM 
AMSTA  SF 

WARREN  MI  48397-5000 

COMMANDER 
WATERVLIET  ARSENAL 
SMCWV  QAE  Q 
B  VANINA 
BLDG  44 

WATERVLIET  NY  12189-4050 

COMMANDER 
WATERVLIET  ARSENAL 
SMCWV  SPM 
T  MCCLOSKEY 
BLDG  253 

WATERVLIET  NY  12189-4050 

TSM  ABRAMS 
ATZKTS 
SJABURG 
W  MEINSHAUSEN 
FT  KNOX  KY  40121 
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10  BENET  LABORATORIES 
AMSTAARCCB 
R  FISCELLA 
G  D  ANDREA 
M  SCAVULO 
G  SPENCER 
P  WHEELER 
K  MINER 
J  VASILAKIS 
G  FRIAR 
RHASENBEIN 
AMSTA  CCB  R 
S  SOPOK 

WATERVLIETNY  12189-4050 

3  ARMOR  SCHOOL 
ATZKTD 
RBAUEN 
JBERG 
A  POMEY 
FT  KNOX  KY  40121 

2  HQ  IOC  TANK 

AMMUNITION  TEAM 
AMSIO  SMT 
R  CRAWFORD 
W  HARRIS 

ROCK  ISLAND  IL  61299-6000 

1  DIRECTOR 
USARMYAMCOM 
SFAE  AV  RAM  TV 
D  CALDWELL 
BLDG  5300 

REDSTONE  ARSENAL  AL 
35898 

2  COMMANDER 

US  ARMY  AMCOM 
AVIATION  APPLIED  TECH  DIR 
J  SCHUCK 

FT  EUSTIS  VA  23604-5577 

1  US  ARMY  CERL 

RLAMPO 

2902  NEWMARK  DR 
CHAMPAIGN  IL  6 1822 


4  DIRECTOR 

US  ARMY  CECOM 
NIGHT  VISION  & 

ELECTRONIC  SENSORS  DIR 

AMSEL  RDNVCM  CCD 

R  ADAMS 

R  MCLEAN 

AYINGST 

AMSEL  RD  NV  VISP 

E  JACOBS 

10221  BURBECK  RD 
FT  BELVOIR  VA  22060-5806 

2  US  ARMY  CORPS  OF  ENGINEERS 

CERDC 
T  LIU 
CEWET 
TTAN 

20  MASS  AVE  NW 
WASHINGTON  DC  203 14 

1  US  ARMY  COLD  REGIONS 

RSCH  &  ENGRNG  LAB 
PDUTTA 
72  LYME  RD 
HANOVER  NH  03755 

1  SYSTEM  MANAGER  ABRAMS 

ATZKTS 
LTCJHNUNN 
BLDG  1002  RM  110 
FT  KNOX  KY  40121 

1  COMMANDANT 

US  ARMY  FIELD  ARTILLERY 
CENTER  AT  FT  SILL 
ATFS  CD 

LTC  BUMGARNER 
FT  SILL  OK  73503-5600 

1  CHIEF 

USAIC 
ATZB  COM 
LTC  TJ  CUMMINGS 
FT  BENNING  GA  31905-5800 

1  NAVAL  AIR  SYSTEMS  CMD 

J  THOMPSON 
48142  SHAW  RD  UNIT  5 
PATUXENT  RIVER  MD  20670 
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8  DIRECTOR 

US  ARMY  NATIONAL  GROUND 
INTELLIGENCE  CTR 
D  LEITER 
M  HOLTUS 
M  WOLFE 
S  MINGLEDORF 
J  GASTON 

W  GSTATTENBAUER 
R  WARNER 
J  CRIDER 

220  SEVENTH  ST  NE 
CHARLOTTESVILLE  VA  22091 

6  US  ARMY  SBCCOM 

SOLDIER  SYSTEMS  CENTER 
BALLISTICS  TEAM 
J  WARD 

MARINE  CORPS  TEAM 

J  MACKIEWICZ 

BUS  AREA  ADVOCACY  TEAM 

W  HASKELL 

SSCNC  WST 

W  NYKVIST 

T  MERRILL 

S  BEAUDOIN 

KANSAS  ST 

NATICK  MA  01760-5019 

9  US  ARMY  RESEARCH  OFC 
A  CROWSON 

J CHANDRA 
H EVERETT 
J  PRATER 
R  SINGLETON 
G  ANDERSON 
D  STEPP 
D  KISEROW 
J  CHANG 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

1  NAVAL  SURFACE  WARFARE  CTR 
DAHLGREN  DIV  CODE  G06 
DAHLGREN  VA  22448 


1  NAVAL  SURFACE  WARFARE  CTR 

TECH  LIBRARY  CODE  323 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

3  NAVAL  RESEARCH  LAB 

I WOLOCK  CODE  6383 
R  BADALIANCE  CODE  6304 
L  GAUSE 

WASHINGTON  DC  20375 

1  NAVAL  SURFACE  WARFARE  CTR 
CRANE  DIVISION 

M  JOHNSON  CODE  20H4 
LOUISVILLE  KY  40214-5245 

2  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 
R  PETERSON  CODE  2020 
M  CRITCHFIELD  CODE  1730 
BETHESDA  MD  20084 

2  NAVAL  SURFACE  WARFARE  CTR 

U  SORATHIA 
C  WILLIAMS  CD  6551 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817 

1  DAVID  TAYLOR  RESEARCH  CTR 
SHIP  STRUCTURES 

&  PROTECTION 
DEPT  CODE  1702 
BETHESDA  MD  20084 

2  DAVID  TAYLOR  RESEARCH  CTR 
R  ROCKWELL 

W  PHYILLAIER 
BETHESDA  MD  20054-5000 

1  OFC  OF  NAVAL  RESEARCH 

D  SIEGEL  CODE  351 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 


t 
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NO.  OF 

COPIES  ORGANIZATION 

8  NAVAL  SURFACE  WARFARE  CTR 
J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
R  D  COOPER  CODE  G32 
J  FRAYSSE  CODE  G33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGREN  VA  22448 

1  NAVAL  SEA  SYSTEMS  CMD 

D  LIESE 

2531  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22242-5 1 60 

1  NAVAL  SURFACE  WARFARE  CTR 

M  LACY  CODE  B02 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

1  OFC  OF  NAVAL  RES 
J  KELLY 

800  NORTH  QUINCEY  ST 
ARLINGTON  VA  22217-5000 

2  NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 

R  CRANE  CODE  2802 
C  WILLIAMS  CODE  6553 
3A  LEGGETT  CER 
BETHESDA  MD  20054-5000 

1  EXPEDITIONARY  WARFARE  DIVN85 

FSHOUP 

2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

1  AFRL  MLBC 

2941  P  ST  RM  136 

WRIGHT  PATTERSON  AFB  OH 

45433-7750 

1  AFRL  MLSS 

R  THOMSON 
2179  12THSTRM  122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 


NO.  OF 

COPIES  ORGANIZATION 

2  AFRL 

F ABRAMS 

J  BROWN 

BLDG  653 

2977  P  ST  STE  6 

WRIGHT  PATTERSON  AFB  OH 

45433-7739 

1  AFRL  MLS  OL 

L  COULTER 
7278  4TH  ST 
BLDG  100  BAY  D 
HILL  AFB  UT  84056-5205 

1  OSD 

JOINT  CCD  TEST  FORCE 
OSD  JCCD 
R  WILLIAMS 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  29180-6199 

1  DEFENSE  NUCLEAR  AGENCY 

INNOVATIVE  CONCEPTS  DIV 
6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 

1  WATERWAYS  EXPERIMENT 
D SCOTT 

3909  HALLS  FERRY  RD  SC  C 
VICKSBURG  MS  39180 

3  DARPA 

M  VANFOSSEN 
SWAX 

L  CHRISTODOULOU 
3701  N  FAIRFAX  DR 
ARLINGTON.  VA  22203-1714 

2  FAA 

TECH  CENTER 
D  OPLINGER  AAR  43 1 
P  SHYPRYKEVICH  AAR  431 
ATLANTIC  CITY  NJ  08405 

2  SERDP  PROGRAM  OFC 

PMP2 

C  PELLERIN 
B  SMITH 

901  N  STUART  ST  STE  303 
ARLINGTON  VA  22203 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  FAA 

MIL  HDBK  17  CHAIR 
L ILCEWICZ 
1601  LIND  AVE  SW 
ANM  115N 
RESTON  VA  98055 

1  US  DEPT  OF  ENERGY 

OFC  OF  ENVIRONMENTAL 
MANAGEMENT 
P  RITZCOVAN 
19901  GERMANTOWN  RD 
GERMANTOWN  MD  20874-1928 

1  DIRECTOR 

LLNL 

F  ADDESSIO  MS  B216 

PO  BOX  1633 

LOS  ALAMOS  NM  87545 

5  DIRECTOR 

LLNL 

R  CHRISTENSEN 
S  DETERESA 
F  MAGNESS 
M  FINGER  MS  313 
M  MURPHY  L  282 
PO  BOX  808 
LIVERMORE  CA  94550 

1  OAK  RIDGE  NATIONAL 

LABORATORY 
RM  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  37831-6195 

1  OAK  RIDGE  NATIONAL 

LABORATORY 
C  EBERLE  MS  8048 
PO  BOX  2009 
OAK  RIDGE  TN  37831 

1  OAK  RIDGE  NATIONAL 

LABORATORY 
CD  WARREN  MS  8039 
PO  BOX  2009 
OAK  RIDGE  TN  37922 


7  NIST 

RPARNAS 
J  DUNKERS 

M  VANLANDINGHAM  MS  862 1 
J  CHIN  MS  8621 
D  HUNSTON  MS  8543 
J  MARTIN  MS  8621 
D  DUTHINH  MS  8611 
100  BUREAU  DR 
GAITHERSBURG  MD  20899 

1  LOCKHEED  MARTIN  MISSILES 

&  FIRE  CONTROL 
R  TAYLOR 

PO  BOX  650003  M  S  WT  93 
DALLAS  TX  75265-0003 

1  HYDROGEOLOGIC  INC 

SERDP  ESTCP  SPT  OFC 
S  WALSH 

1 155  HERNDON  PKWY  STE  900 
HERNDON  VA  20170 

3  DIRECTOR 

SANDIA  NATIONAL  LABS 
APPLIED  MECHANICS  DEPT 
DIV  8241 
J  HANDROCK 
Y  R  KAN 
J  LAUFFER 
PO  BOX  969 

LIVERMORE  CA  94550-0096 

3  NASA  LANGLEY  RSCH  CTR 

AMSRL  VS 
WELBERMS  266 
F  BARTLETT  JR  MS  266 
G  FARLEY  MS  266 
HAMPTON  VA  23681-0001 

1  NASA  LANGLEY  RSCH  CTR 

T  GATES  MS  188E 
HAMPTON  VA  23661-3400 

1  USDOT  FEDERAL  RAILRD 

M  FATEH  RDV31 
WASHINGTON  DC  20590 
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1 

1 

1 
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FHWA 
E  MUNLEY 

6300  GEORGETOWN  PIKE 
MCLEAN  VA  22101 


COMPOSITE  MATERIALS  INC 
C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 


CENTRAL  INTLLGNC  AGNCY 
OTI  WDAG  GT 
W  L  WALTMAN 
PO  BOX  1925 
WASHINGTON  DC  20505 


2  COMPOSIX 
D  BLAKE 
L  DIXON 
120  0  NEILL  DR 
HEBRUN  OHIO  43025 


MARINE  CORPS  INTLLGNC  ACTVTY 
D  KOSITZKE 

3300  RUSSELL  RD  STE  250 
QUANTICO  VA  22134-501 1 

DIRECTOR 

NATIONAL  GRND  INTLLGNC  CTR 
LANG  TMT 
220  SEVENTH  ST  NE 
CHARLOTTESVILLE  VA 
22902-5396 

DIRECTOR 

DEFENSE  INTLLGNC  AGNCY 
TA  5 

K  CRELLING 
WASHINGTON  DC  20310 

GRAPHITE  MASTERS  INC 
J  WILLIS 

3815  MEDFORD  ST 

LOS  ANGELES  CA  90063-1900 

ADVANCED  GLASS  FIBER  YARNS 
T  COLLINS 

281  SPRING  RUN  LANE  STE  A 
DOWNINGTON  PA  19335 

COMPOSITE  MATERIALS  INC 
D  SHORTT 
19105  63  AVENE 
PO  BOX  25 

ARLINGTON  WA  98223 

COMPOSITE  MATERIALS  INC 
R  HOLLAND 
11  JEWEL  CT 
ORINDA  CA  94563 


4  CYTEC  FIBERITE 

R DUNNE 
D  KOHLI 
M  GILLIO 
RMAYHEW 
1300  REVOLUTION  ST 
HAVRE  DE  GRACE  MD  21078 

2  SIMULA 

J  COLTMAN 
RHUYETT 
10016  S51ST  ST 
PHOENIX  A Z  85044 

1  SIOUX  MFG 
BKRIEL 
PO  BOX  400 

FT  TOTTEN  ND  58335 

2  PROTECTION  MATERIALS  INC 
M  MILLER 

F  CRILLEY 
14000  NW  58  CT 
MIAMI  LAKES  FL  33014 

3  FOSTER  MILLER 
J  J  GASSNER 

M  ROYLANCE 
WZUKAS 
195  BEAR  HILL  RD 
WALTHAM  MA  02354-1196 

1  ROM  DEVELOPMENT  CORP 

ROMEARA 
136  SWINEBURNE  ROW 
BRICK  MARKET  PLACE 
NEWPORT  RI 02840 
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NO.  OF 

COPIES  ORGANIZATION 

2  TEXTRON  SYSTEMS 

T  FOLTZ 
M  TREASURE 
201  LOWELL  ST 
WILMINGTON  MA  08870-2941 

1  JPS  GLASS 

L  CARTER 
PO  BOX  260 
SLATER  RD 
SLATER  SC  29683 

1  O  GARA  HESS  &  EISENHARDT 
M  GILLESPIE 

9113  LESAINT  DR 
FAIRFIELD  OH  45014 

2  MILLIKEN  RSCH  CORP 
H  KUHN 

M  MACLEOD 
PO  BOX  1926 
SPARTANBURG  SC  29303 

1  CONNEAUGHT  INDUSTRIES  INC 
J  SANTOS 

PO  BOX  1425 
COVENTRY  RI 02816 

2  BATTELLE  NATICK  OPNS 
J  CONNORS 
BHALPIN 

209  W  CENTRAL  ST  STE  302 
NATICK  MA  01760 

1  BATTELLE  NW  DOE  PNNL 

T  HALL  MS  K231 
BATTELLE  BLVD 
RICHLAND  WA  99352 

1  ARMTEC  DEFENSE  PRODUCTS 
S  DYER 

85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 

2  GLCCINC 
J  RAY 

M  BRADLEY 

103  TRADE  ZONE  DR  STE  26C 
WEST  COLUMBIA  SC  29170 


NO.  OF 

COPIES  ORGANIZATION 

3  PACIFIC  NORTHWEST  LAB 

M  SMITH 
G  VAN  ARSDALE 
R  SHIPPELL 
PO  BOX  999 
RICHLAND  WA  99352 

2  AMOCO  PERFORMANCE 
PRODUCTS 

M  MICHNO  JR 
J  BANISAUKAS 
4500  MCGINNIS  FERRY  RD 
ALPHARETTA  GA  30202-3944 

1  SAIC 

M  PALMER 

2109  AIR  PARK  RD  S  E 
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